Dislocation-mediated growth of bacterial cell walls 
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Recent experiments have illuminated a remarkable growth mechanism of rod-shaped bacteria: 
proteins associated with cell wall extension move at constant velocity in circles oriented approxi- 
mately along the cell circumference (Garner et al., Science (2011), Dominguez- Escobar et al. Science 
(2011), van Teeffelen et al. PNAS (2011). We view these as dislocations in the partially ordered 
peptidoglycan structure, activated by glycan strand extension machinery, and study theoretically 
the dynamics of these interacting defects on the surface of a cylinder. Generation and motion of 
these interacting defects lead to surprising effects arising from the cylindrical geometry, with impor- 
tant implications for growth. We also discuss how long range elastic interactions and turgor pressure 
affect the dynamics of the fraction of actively moving dislocations in the bacterial cell wall. 
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Bacterial cell walls are composed of peptidoglycan 
(also called murein) , which endows them with shape and 
rigidity. Their architecture and growth have been the 
subject of active research for many decades, in particular 
for gram-negative bacteria whose cell walls consist of sin- 
gle or few layers of glycan strands crosslinked by peptides 
0-13: 0- While some models assume a very ordered struc- 
ture, recent experimental work Q suggests the structure 
is more disordered. We view the peptidoglycan mesh as 
a partially ordered two-dimensional crystal with a large 
number of defects to account for the disorder in the struc- 
ture. The rod shape of many bacteria (e.g.: Escherichia 
coll), together with mutant variants that grow but fail 
to complete cell division, make a cylindrical geometry a 
natural one to study. To easily add material to this or- 
dered structure, one must clearly break the periodicity 
and create a defect in the structure. An especially im- 
portant class of defects are termed dislocations, known 
to be important in determining the mechanical properties 
of metals and other crystalline or polycrystalline materi- 
als, such as their strength and plasticity y . Dislocations 
are known to have long-ranged elastic interactions, which 
are logarithmic in the distance (but not isotropic, as is 
the case for vortices in superfluids). The "elementary 
charge" of these topological defects is the Burgers vector 
6, which is often at right angles to the direction of inser- 
tion of the new strand of material, and is a lattice vector 
of the structure. Here, we apply dislocation theory to 
the problem of bacterial growth, proposing a simplified 
model inspired by the elongation of bacterial cell walls, 
which we are able to solve both analytically and via com- 
puter simulations. Fig. 1 illustrates an idealized picture 
of defects in a cylindrical geometry. For simplicity, we 
show a square lattice with lattice vectors parallel and 
perpendicular to the cylinder's long axis., although the 
actual peptidoglycan mesh is rectangular (and in fact the 
lattice vectors might have a non-trivial angle with respect 
to the cylinder's axis 0)- 

The model treated below will in fact consist of a large 
number of dislocations, and so the structure we treat is 



far from a perfect crystal. Working with defects of a 
crystal provides a convenient and numerically efficient 
method to take the disorder into account. We expect 
that the unit cell of Fig. 1 typically contains two gly- 
can strands, as it is only this larger unit cell that re- 
spects the local crystalline symmetry; see for example 
Refs. [1, 0. Recent experiments on both gram- negative 
[lo| and gram-positive bacteria 11, 12 1 track fluores- 
cently labeled proteins known to correlate strongly with 
the addition of peptidoglycan subunits, and have shown 
that these proteins move at roughly constant velocity, 
approximately along the cylinder's circumference. We 
view these strand extension centers as edge dislocations 
in the ordered structure, with a Burgers vector oriented 
along the cylinder's long axis (the direction of the Burg- 
ers vector depends on the direction of insertion of the 
new strand, see Fig. 1). Extending the end of an in- 
serted strand {i.e., the core of an edge dislocation) in- 
volves breaking peptide bonds to allow extra sugar units 
into the lattice, together with additional short peptide 
crosslinks 0] . In dislocation theory, this type of motion 
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FIG. 1: Schematic illustration of active (arrows) and inactive 
(asterisk) dislocations in an otherwise ordered peptidoglycan 
mesh. The dislocations with arrows attached are activated 
by the enzymatic machinery and move with velocity v. Those 
with asterisks are inactive. 
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is referred to as dislocation climb. In the following, we 
treat the protein motion and the dislocation motion syn- 
onymously, assuming that the motion of the MreB and its 
associated enzymes is fully correlated with the insertion 
of new material into the cell wall. This idea was intro- 
duced in 0, and here we develop it further and deduce 
various biological insights and predictions. As pointed 
out by Burman and Park [isj , glycan strand extension is 
somewhat analogous to the action of a DNA polymerase. 
Our dislocation perspective allows us to take into account 
long range elastic interactions between the "murein ex- 
tension centers" of Ref . [l^ . 

Although extensive work has been done on dislocation 
theory over the last century this biophysics prob- 
lem is quite different from materials science and con- 
densed matter physics in several respects: the climb of 
dislocations necessarily involves the exchange of mate- 
rial, which in 3d crystals involves diffusion of interstitials 
or vacancies. Hence, except at high ternperatures, dis- 
location glide dominates the dynamics [iJ]. Here, dis- 
location climb is the central process mediating cell wall 
growth, with sugars and amino acids essential for the 
climb process synthesized and arriving from the interior 
of the bacteria. A second, obvious, difference, regards 
the unusual cylindrical geometry which we study here. 
This feature leads to a number of interesting properties, 
such as the exponential decay of dislocation interactions 
along the cylinder's long discussed below. Dis- 

locations in a cylindrical geometry were considered in a 
very different biophysics problem in Ref. [isj . studying 
tail-sheath contraction in a bacteriophage. A final dif- 
ference from conventional materials science is that here 
the dislocation climb itself alters the geometry, as it is 
this very motion that grows the cell wall along the axis 
of the cylinder. This idea can be used to estimate the 
number of actively climbing dislocations per bacterium, 
using experimental data: The typical time of division of 
the bacteria is of the order of tens of minutes, in which 
the bacteria elongate by approximately 1 ^m. The lat- 
tice spacing of the peptidoglycan along the long axiSj6, 
is believed to be of the order of several nanometers [7| . 
If we assume the area of the hemispherical end caps re- 
mains fixed, the dynamics of the cylindrical surface area 
S(t) (taking all dislocations as completely independent, 
for now), is given by: 



tive labeling. Current technology does not yet enable a 
direct determination of Nac, since only a subset of the to- 
tal number of active dislocations is fiuorescently marked 

iMli. 



dS_ 

dt 



= NacVb, 



(1) 



where Nac is the total number of actively moving dis- 
locations and V is their velocity. Taking the measured 
velocity v to be several tens of nanometers/second (lol - 
12l | , we find that a few tens of active dislocations moving 
along the circumference and growing the cell wall would 
be sufficient to account for the measured growth rate. 
This estimate is consistent with pioneering work of Ref. 
[isj . obtained using a very different method of radioac- 



The structure of the manuscript is as follows. We first 
define the model and the approximations made, and de- 
scribe the force on a dislocation, including the elastic 
interactions between dislocations on a cylinder. The im- 
plications of the large turgor pressure inside the bacteria 
and the effects of the long ranged elastic interactions are 
highlighted. We then propose a novel set of equations 
for the growth dynamics and its coupling to the num- 
bers of active and inactive dislocations, and show that 
they often lead to exponential lengthening of a single 
bacterium. The exponential growth rate itself depends 
on a few simple parameters, with a well defined micro- 
scopic interpretation. The theoretical expectations are 
illustrated and visualized by numerical simulations, both 
in the main text and the Supplementary Material. A nu- 
merical calculation presented in the Supplementary Ma- 
terial is used to estimate the disorder strength due to 
the elastic interactions with a large number of randomly 
positioned dislocations, chosen according to the biologi- 
cal parameters. Certain parameters of the rate-equations 
model (Eqs. ([5]) and ([6])), can not be determined from 
a numerical simulations, and are dictated by the under- 
lying biochemistry, for example v(G) (the dependence of 
the dislocation velocity on the driving force) and the dis- 
location processivity 71. 



THE MODEL 

Interesting recent works have modeled and simulated 
the peptidoglycan structure in molecular detail, and have 
predicted a number of non-trivial phenomena which were 
indeed observed experimentally fwi Kl\ . Our approach is 
somewhat different (and perhaps complementary), since 
we study a highly simplified model and focus on a rel- 
atively dilute concentration of defects in an otherwise 
regular peptidoglycan mesh. As a result both analytic 
and numerical progress is possible, although we suppress 
fine details of the structure. 

Consider, motivated by bacterial cell wall growth, a 
number of interacting edge dislocations with a Burgers 
vector b oriented in both directions along the cylinder's 
long axis, as illustrated schematically in Fig. 1. The ra- 
dius of the cylinder, R, will be assumed to be constant, as 
is approximately true for rod-shaped bacteria. For sim- 
plicity, consider an infinite cylinder, neglecting the effects 
of the boundaries (Experiments indicate that strand ex- 
tension dynamics does not change significantly near the 
cylinder's end caps, consistent with this assumption ^18:]). 

We also simplify by assuming linear elasticity, except 
near dislocation cores. Various works have indicated, 
both from the theoretical and experimental per- 
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spcctivc, that the high osmotic pressure [s, 21|, referred 
to as turgor pressure in bacteria and plant cells, pro- 
duces large strains compared to the relaxed state. We 
shall take the elastic moduli as constants, which in prac- 
tice should be considered as effective constants relative 
to the working turgor pressure. 

In order for a dislocation to be able to "climb" (i.e., 
move perpendicular to its Burgers vector ^ and con- 
tribute to the growth process), various proteins have to 
be present at the dislocation core, responsible for "re- 
cruiting" and assembling the sugars and peptides nec- 
essary for the construction of the peptidoglycan mesh. 
Hence, we separate the dislocations into two populations, 
active and inactive ones. Inactive defects still exert elas- 
tic forces on the active ones, and create an effective dis- 
ordered energy landscape for them. They also represent 
favorable locations for the creation of new active dislo- 
cations: the elongation machinery can attach to the free 
strands (i.e., the inactive dislocations), creating active 
ones. The finite processivity of the elongation machinery 
will also give rise to the opposite process, whereby an 
active dislocation can become inactive when the machin- 
ery falls off. In general, the rates of these processes will 
not balance, since the system is never in steady-state but 
is instead constantly growing via dislocation climb, as it 
incorporates material from the third dimensions, a situa- 
tion rarely encountered in conventional materials science 
i- 



The force on a dislocation 

Elastic stresses exert a Peach-Koehler force 0, Q on 
a dislocation (analogous to the Magnus force acting on 
vortices which in our case is given by [il ): 



F„ 



(2) 



where a is the two-dimensional stress tensor of the pep- 
tidoglycan mesh, and we assume a Burgers vector along 
the X-axis. 

For a cylinder of radius R we have ciyy = 2axx ~ pR , 
where x denotes the coordinate along the long axis, and 
p is the turgor pressure, and Fy = —^bpR. In addition to 
the contribution of the turgor pressure, if the free energy 
is changed by U by the biochemical process of adding one 
unit cell to the peptidglycan mesh, it will contribute an 
additional force of U/b in the y direction, as can be seen 
using the principle of virtual work. 

Under physiological conditions, we expect that the dis- 
locations are in the overdamped regime, with the dis- 
location velocity proportional to the force, i.e., Vi = 
'^jfJ'ijFj, where fiij is a mobility tensor with glide and 
climb components that depends on the extension machin- 
ery and the abundance of sugars, peptides, etc. In the 



following we assume that fiij is diagonal, with fx^x and 
fiyy describing glide and climb mobilities, respectively. 
We expect that the mobility tensor itself will have a tur- 
gor pressure dependence, as the resulting forces can lower 
the activation barriers of the various biochemical path- 
ways involved in the process. Thus, the observed velocity 
of the strand extension machinery should depend on the 
excess pressure for rod-shaped bacteria. 



The role of interactions 

In condensed matter physics, the long-range elastic in- 
teractions between dislocations can have important con- 
sequences, and have recently been suggested to lead to 
glassy effects and non-thermal, heavy-tailed, dislocation 
velocity distributions [251427 1. 

To illuminate their importance in the biological con- 
text, we have solved for the interaction energy of two 
dislocations with antiparallel Burgers vectors zLb on an 
infinite cylinder, separated by a distance x along the long- 
axis and y along the circumference [28| . 



^, , Ab'^h. X / y \ . 1 / X , , rX 

Fix.y) = \i — CSC — Ismh — Icsch — 

^ 2 ^ 2R ^2W ^2R' ^ 2R 

2R X - iv 
+log(— sinh[^]]+C.C. 



(3) 



Here, R the radius of the cylinder, h the cell wall thick- 



ness, and A 



27r(i-i/) 



, with /i the shear modulus and 



V the Poisson ratio (we assume isotropic elasticity the- 
ory, for simplicity). For distances x ^ R along the axis 
of symmetry, Eq. (3)) shows that the interactions fall 
off exponentially [29[. For distances much smaller than 
the cylinder radius, the interaction energy reduces to its 
form in two-dimensions, as expected. Clearly, the rele- 
vant scale for the interaction energy is Ab^h. Fig. la 
in the Supplementary Materials shows an example of the 
equal energy contours. Close to the origin, taking a cut 
parallel to the x axis would give a graph with two min- 
ima, corresponding to the two dislocations sitting at a 45 
degree angle - this is expected, since the physics close to 
the origin should be insensitive to the finite cylinder ra- 
dius: this arrangement is indeed the stable configuration 
of two dislocations in Euclidean space, when climb pro- 
cesses arc prohibited [30t . Fig. lb in the Supplementary 
Materials illustrates how interactions affect an activated 
dislocation (i.e., one with strand-extending machinery 
attached) attempting to move upwards from a fixed, un- 
activated defect. 

Let G represent the combined effect of the chemical 
forces provided by enzyme-mediated polymerization and 
those caused by turgor pressure, as described by Eq. 
(I2|). The sign of this force depends on the sign of the 
Burgers vector, and acts in the direction which corre- 
sponds to elongation of the existing strands. For a ver- 
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tical separation y in the range b < y <^ 27ri?, we can 
use Eq. ([3]) and add a term Efieu ~ ^Gy correspond- 
ing to the force G. This leads to the potential energy 
Uiy) = 2E^ + (Ab'^h) log(y/6) - Gy, where 2E^ is the 
energy at separation y = b and G = 0. The logarith- 
mic interaction term clearly leads to a long range elastic 
force that falls off slowly, like 1/y. Upward motion will 
require thermal activation unless y > y* = Ab^h/G. Es- 
cape from this logarithmic binding will be relatively easy 
for G large enough to insure y* < b, i.e., G > G* = Abh . 
Healthy, growing cells in rich nutrient environments are 
likely to be in this regime and thus escape the confining 
effect of clastic forces. Elastic interactions will neverthe- 
less become important when G is diminished by osmotic 
shock, nutrient deprivation, or by application of antibi- 
otics that interfere with the powerful strand extension 
machinery. 

Fig. 2 illustrates the dynamics of several tens of ac- 
tive dislocations under the influence of a circumferential 
"field" or driving force G on each dislocation with pair- 
wise elastic interactions given by Eq. ([3|), moving in the 
disordered background created by inactive dislocations, 
to be elaborated upon below. For a weakly disordered 
system, one obtains ordered circumferential motion (see 
Fig. 3 in the Supplementary Material) , while for a disor- 
der level chosen according to the biological parameters, 
the elastic interactions are strong enough to occasionally 
pin the motion. This could be related to pauses in the cir- 
cumferential motion observed in recent experiments . 
The condition for pinning to be important follows from 
comparing the typical elastic interaction force with the 
nearest dislocation a distance r„„ away (which is highly 
likely to be inactive) and the net driving force G. Pinning 
is expected to occur when: 

^ = Ahb^n'/^ > G, (4) 

where n is the dislocation arcal density. With parameter 
estimates for E. coli (see the section "Estimation of the 
biological parameters'" in the Supplementary Material) 
we find that the ratio x = Ahb^n^/^ /G is of the order of 
unity, suggesting that pinning events are likely to occur. 
However, due to the finite processivity of the active dislo- 
cations, pinning can be bypassed by the growth dynam- 
ics, as the elongation machinery can drop off a pinned 
dislocation and continue the strand elongation elsewhere. 
However for biological conditions (nutrient depiction or 
low turgor pressure) such that % ^ 1, pinning events 
could become much more important. 

Evidently, experimental measurements of interaction- 
induced correlations between moving dislocations could 
provide a useful probe of the elastic properties of the 
peptidoglcan. Estimates of the relative strength of the 
interactions compared to the other forces involved in the 
problem are provided in the Supplementary Material, 
and discussed below. 



DISLOCATION DYNAMICS 

Eq. ([T]) expresses the rate of change of the bacterial 
surface area S{t). To proceed further, we must account 
for the dynamics of the number of active and inactive dis- 
locations, controlled by a variety of interesting biological 
processes: 

1. Dislocation pair creation: it is plausible that the 
nucleation of new growth sites should be proportional to 
the cylinder's area (see Fig. 3b). For growth under rich 
nutrient conditions, we expect the area density of pcp- 
tidoglycan growth sites {i.e., activated dislocations) to 
be constant, which requires a source of new dislocations 
as the bacteria elongates with fixed radius. Dislocation 
pairs could in principle be created by thermal activation 
1131, although this process is likely to be rare in a 
covalently bonded peptidoglycan mesh. Alternatively, a 
dislocation pair can form by breaking bonds in the pepti- 
doglycan mesh, via an enzymatic process or via external 
damage by radiation or severe stresses [slj . 




0.5 3 



FIG. 2: Numerical simulation of the dynamics of 30 active 
dislocations, with simulation parameters matching the bio- 
logical ones. The dislocations are driven by a combination 
of chemical forces and turgor pressure in the circumferen- 
tial direction, proportional to the sign of the Burgers vec- 
tor. 10,000 inactive dislocations create a disordered potential 
for the motion of the active dislocations. The temperature 
is zero and ^^^^ — 450. The x and y coordinates of the 
dislocations were chosen randomly and uniformly along the 
axes, and W/ L = 2lS — 1. The climb to glide mobility ratio 
is fic/f^g ~ 10. The red and blue points mark the starting 
positions of the ±6 dislocations, and the red and blue lines 
correspond to their trajectories. The black circles mark the 
end of the cylinder. A snapshot of the simulation is taken 
after t = 0.7^. The numerous inactive dislocation (marked 
as black dots) create a disordered energy landscape which can 
trap the active dislocations and inhibit further growth. In the 
Supplementary Material, further details regarding the simula- 
tion are given. 
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2. The finite processivity of eiongation enzymes: th.e 
elongation macliinery can fall off an active dislocation, 
and render it inactive. In an elegant recent work 33[, 
finite processivity was implicated in straightening of an 
initially curved Caulobacter crescentus, where the work- 
ing hypothesis was that the growth occurs via partial 
"hoops" extending perpendicular to the lon g ax is of the 
rod shaped bacteria, as described in Refs. [10|-ll2|. For 



C. crescentus, processivity was such that on average 1/5 
of the circumference is traversed before the machinery 
falls off [3^ . It would be interesting to perform a similar 
study for E. coli or Bacillus subtilis. Growing filamentous 
E. coli in curved microchambcrs was also demonstrated 
[s^ l- proving the feasibility of such an experiment for E. 
coli. 

3. The activation of an inactive dislocation; the re- 
verse process assembles the elongation machinery onto 
an existing (static) glycan strand end and turns it into 
an actively growing dislocation. The rate of this process 
is proportional to the number of inactive dislocations. 

4. Annihilation processes: two dislocations, both ac- 
tive or one inactive, can "collide" , completing a full cir- 
cumferential strand. The elongation machinery could 
then immediately continue to move and grow the cell wall 
(as in Fig. 3a), in which case the annihilation process 
plays no role. Alternatively, it could tie the strands to- 
gether, dislodge and become available for a new inactive 
dislocation. In complete analogy to the theory of chemi- 
cal reactions, the rate of this latter process per unit area 
is proportional to the product of the density of the active 
dislocation and that of the total number of dislocations. 



Upon collecting together these four processes, we de- 
scribe the dynamics of the numbers of active and inactive 
dislocations on a cylindrical wall of surface area S(t) by 
supplementing Eq. ([1]) with the following: 



dNg, 
dt 



dNi, 



= TS - 71 iVac + llNin - IzvNaciriac + Tlin) , (5) 



= r5 + "flNac - "f2Nin + 2745 - JsVriacNin. (6) 



Nac and Nin denote the number of ac- 

^ Nac/S and 



dt 

In the above 

five and inactive dislocations, and 
nin = Nin/ S denote their densities. The term F arises 
from the nucleation of dislocation pairs, one active and 
one inactive, by a process like that shown in Fig. 3b, 
while 7i represents the processivity and 72 describes the 
activation of an inactive dislocation (see Fig. 3a). The 
rate 73 represents the annihilation of a dislocation pair 
(which is typically proportional to the velocity Q). Fi- 
nally, 74 describes creation of tightly bound inactive dis- 
location pairs by an enzyme which adds to the reservoir of 
activatable dislocations by cutting covalent bonds (Fig. 
3a) . Note that no spatial structure is present in the above 
equations: the dynamics is "mean- field" like, i.e., it as- 
sumes that the dislocations are well-mixed and that their 
density is approximately uniform on the cylinder [ssj . A 
more sophisticated theory would explore the effects of 
spatial heterogeneities on the dynamics [Hj. 

We search for steady state densities Uac and n.i„ with 
these equations, taking into account Eq. ([T]). The steady 
state dislocation densities then obey: 



(a) 
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(b) 



FIG. 3: (a) Active dislocation pulled past an inactive one by 
the enzymatic machinery. In this case an initial (inactive) 
dislocation pair could be created by an enzyme that cuts gly- 
can bonds (rate constant 74). The elongation machinery then 
assembles around an inactive dislocation turning it into an 
active one (rate constant 72). (b) Direct insertion of a glycan 
strand fragment to create an active and inactive dislocation 
pair (rate constant F). 



rii^vb = r - 7inac + 72'^m - ld.vnac{nac + riin), (7) 



ninUacVb = r + 7i71qc - 72?^^n + 274 - Jsvriacnin- (8) 

The terms on the left arises because the dislocations 
are automatically diluted as the cell wall grows. From 
Eq. ([T]), we see that the length of the bacteria L{t) ~ 
grows exponentially in the steady state. 



vbt 



(9) 



While we cannot rule out the possibility of cell cy- 
cle effects in the elongation rate, our simplified model is 
consistent with the existing data for exponential growth, 
for bacteria growin g in one dimensional channels with 
abundant nutrients [36[. As discussed below, the G force 
driving dislocation climb enters our dynamical equations, 
through (1) F(G), the escape rate of bound dislocations, 
(2) the velocity v = v{G), (3) through the capture cross- 
section embodied in 73 = 73(G), and (4) the creation 
rate of new dislocations. In particular F and 74, might 
depend on G, as it could change the barriers for thermal 
activation (see 37[ for a related discussion). 
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THE BIOLOGICAL PARAMETERS 

In the Supplementary Material, we estimate the rele- 
vant biological parameters for the case of E. coli, in or- 
der to elucidate which of the previous physical processes 
should be important. We estimate the energy of cre- 
ating a new dislocation pair using the measured elastic 
constants, and find that it is much larger than the tem- 
perature. Using the measured turgor pressure, we find 
that the force exerted on each of the dislocation by the 
stresses arising from it is similar in strength to those due 
to the biochemical reaction involved in the elongation 
process. This suggests that changing the turgor pressure 
or alternatively, putting the bacteria under mechanical 
stresses, can result in significant changes in the velocity 
of the dislocations. The density of inactive dislocations 
can be estimated from the strand length distribution, and 
is found to be much larger than that of active disloca- 
tions. Finally, the typical elastic interaction forces be- 
tween nearby active dislocations are negligible compared 
to the driving force in the circumferential direction, but 
the elastic interactions between the active dislocations 
and the numerous inactive ones are comparable to the 
driving force, suggesting that the possibility of disloca- 
tion pinning, discussed earlier, could occur. 

BACTERIAL ELONGATION RATE 

The assumption ni„ ^ Uac simplifies Eqs. [5] and |6l 
Wc explore here two possible scenarios: 

1. A two-stage process: An enzyme cuts glycan links 
in the mesh to produce a pair of inactive dislocations, 
while other proteins recruited by the MreB then acti- 
vate a small subset of them, accounting for the small 
subset of active dislocations (see Fig. 3a). Thus far, 
the annihilation of dislocations does not seem to be ob- 
served in experiments [isj , suggesting that we could also 
ignore 73. With these assumptions, we can simplify Eqs. 
([5|) and ([6]) in the steady state: n^^ufo = — linac, 

ninHacVb = JlUac - + 274. 

Upon adding these equations and using riin 3> n-ac we 
obtain: ninUacvh = 274 , which when combined with the 
previous equations leads to a cubic equation for tIqc, the 
density of active dislocations: 

riac&V+7inL6« = 27274. (10) 

For 6^7274 S> 7i (short processivity of the elongation 
machinery) we find Uac = (^5^)^, while for &W7274 ^ 
ll get n,, = 

It is useful to rewrite these scaling relations in terms 
of the exponential growth rate A = Uacvh, where L{t) = 
Lqc^* . This leads to the relations: A ~ {2bv^274)^ , 



for short processivity, and: A = y ^^^^'^ in the oppo- 
site regime. An experimentally accessible control {e.g.: 
temperature) for which one of the microscopic processes 
(e.g.: 71) is much more sensitive than the others, would 
provide a way to test these scaling relations without ex- 
perimentally measuring all of the microscopic rates. 

2. Direct insertion of a nascent strand: Another pos- 
sible scenario is that dislocation pairs are formed by in- 
serting material into the mesh (but on new sites, between 
the existing mesh, see Fig. 3b), which are further acted 
upon by the elongation machinery to produce the growth. 
If this is the dominant mechanism we can neglect 74. 
If in addition we neglect the annihilation term, as dis- 
cussed previously, we obtain the steady state equations: 
nl^vb = r-7inac + 72n„, n^nTlacVb = r + 7inac-72"-m- 

Upon again invoking ^ 1, wc find that nacninbv 
2r, which leads to a quadratic equation for Uac- 

7i"L - T'^ac - = 0. (11) 

For 71 » we obtain n,, « 

Another useful limit can be obtained by setting 72 = 0. 
In this limit active dislocations can only be created by 
the direct insertion, and inactive dislocations to not play 
a role. In this case, which was studied in Ref. the two 
equations decouple, and one obtains a quadratic equation 

for Hac- 

We anticipate that experiments will soon determine 
the density of active dislocations, thus allowing tests of 
the above scaling relations. As explained previously, pa- 
rameters such as the processivity could be experimentally 
determined independently, allowing additional checks. 
Another interesting biological scenario is one where the 
creation of new dislocations is abruptly stopped, by dis- 
abling the relevant protein. In this case 74 = F = 0, 
and subsequent steady-state growth is linear rather than 
exponential, since the total number of dislocations is con- 
stant in this case. From Eqs. (HI and ^ we can deter- 
mine the typical transient time r for the growth to cross 
over to the linear regime is, r = — } — . This and re- 

o ' 7l-t-72 

lated experiments could provide further insight into the 
microscopic rate constants discussed above. 

Protein regulation and exponential growth.- Cell wall 
growth by itself is unlikely to be a bottleneck process that 
limits bacterial growth [38[. In fact, bacterial cells prob- 
ably regulate the concentration of all proteins, DNAs, 
RNAs etc. involved in the cell growth, to yield a coordi- 
nated set of doubling times determined by richness of the 
media, temperature etc. In our model, parameters such 
as the dislocation velocity, driving force G and density 
of active dislocations will be determined by the regula- 
tory network that controls cell growth. One important 
feedback mechanism is surely the "geometrical dilution" : 
as the cell wall grows, the concentrations of all the pro- 
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tcins in the cells are diminished. If we assume ribosome 
production is the limiting factor of the growth jsoj (ribo- 
somes constitute of most of the dry weight of the cell), 
and the relevant proteins are maintained at constant con- 
centration, the cell wall dynamics discussed here would 
be consistent with the ribosomal constraints. Regula- 
tion of concentration is appealing, as it is a local and 
robust mechanism, consistent with the biological ideas 
about self-organization Locality is also connected 

with the appearance of the surface area in Eqs. ([Sj and 
([5]). as discussed earlier. Even if some more complicated 
biological process is involved, we expect predictions such 
as Eqs. ((TT|) to hold; one should be able to check them 
directly by measuring the rates of the underlying indi- 
vidual processes, and by engineering sudden changes in 
rates such as 71, 72 or 74. Finally, we note that the 
aforementioned geometric dilution is quite generic, and 
will appear in every case where the cell volume grows. 
In Ref. [4l|, section 2.4, this dilution is accounted for 
by adding a term to the protein degradation. In our 
case, however, the dilution of the active dislocations is 
not proportional to their concentration (as is usually the 
case with protein degradation), but enters quadratically. 
This comes about since the active dislocations themselves 
are responsible for the change in surface area. 



CONCLUSIONS AND OUTLOOK 

In this paper we used the continuum elastic theory to 
account for recent experiments on the growth of single 
bacteria. Our theoretical framework exploits the par- 
tially ordered structure peptidoglycan mesh to obtain an 
economical description of the growth dynamics in terms 
of dislocation defects. The forces acting on these de- 
fects lead to predominantly climb dynamics, quite dif- 
ferent from those typically found in condensed matter 
physics and materials science. Our results suggest that 
osmotic pressure could affect the measured dislocation 
velocities, and that the elastic interactions could lead to 
pinning of the dislocation motion. Consideration of var- 
ious processes associated with dislocation motion, and 
their creation and annihilation, leads to a coupled set of 
ordinary differential equations that govern the dynamics 
of the number of active and inactive dislocations, and the 
overall elongation rate. Solving these under conditions of 
abundant nutrients leads to exponential growth, consis- 
tent with the experimental observations. Transients and 
other modes of growth are also discussed. In the future, 
it would be interesting to compare and contrast bacterial 
elongation with pollen tube growth, another case where 
the growth occurs on a cylinder with an approximately 
constant radius 



cells. Another natural extension of this approach might 
be to model the growth of approximately spherical cells 
, such mutant E. coli. Finally, it would be useful to ex- 
tend this theory to the thicker cell walls of gram-positive 
bacteria, which experimentally show similar phenomenol- 

ogy ^m- 
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42l |. One might also speculate that pla- 



nar arrangements of cells in, say, embryos, grow by cell 
divisions sited preferentially in dislocation cores. The 
resulting dislocation climb takes place in array of living 
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Supplementary Material 



Dislocation interactions on a cylinder 

Consider edge dislocations on the surface of a cylinder, 
with radius R and infinite length. Let us denote the 
coordinate along the symmetry axis by x, and the other 
by y, i.e., y is periodic with a period of = 2ttR. We 
shall assume throughout that the Burgers vector of all 
dislocations point in the ±x direction, with magnitude b. 

We will not consider out-of-plane buckling, which is 
suppressed by the high turgor pressure inside the bac- 
teria [H. Confined to the surface of the cylinder, the 
Laplacian operator is equivalent to that in an infinite 
two-dimensional space together with the periodicity re- 
quirement. Therefore, to find the stress exerted by a dis- 
location at position (xq, i/q) on another at position {x, y), 
we have to solve the bi-harmonic equation for the Airy 
stress function Q on a strip with periodic boundary con- 
ditions. 

The stresses without the periodicity requirement are 
given by (up to a sign depending on the direction of the 
Burgers vector): 



The force on an additional dislocation at (x, y) due to 
this stress will then be given by the Peach-Koehler force 



P _ 1 cylinder 



(7) 



(8) 



Although the stress component <jyy does not enter the 
forces in this geometry, we will also evaluate it, for com- 
pleteness. 

The sums of Eqs. ^ can be evaluated using 
the Sommerfeld- Watson transformation We illus- 

trate this method for the first sum. Since the function 
g{z) = cot(7rz) has only simple poles of residue unity that 
lie on the x axis at integer values, the sum in Eq. ([1]) can 
be written as the complex contour integral: 



g{z)f{z)dz, 



(9) 



with: 



- ~Abh{y - y^mx - xof + [y - yof]/r\ (1) 



Ahh{y - yo)[{x - xof - {y - y^f]lr\ (2) 



and 



'xy 



= Abh{x - xq)[{x - xof - {y- yo)^]/f'* 



(3) 



withr^ = {x~xo)'^ + {y—yQ)'^ and A = 27r(i-g-) ' ''^i^h /j, the 
shear modulus and cr the Poisson ratio. Note that these 
stresses are in two-dimensions, and hence their units are 
force per unit length. 

To describe stresses on a cylinder, the dislocation at 
(xQ^yo) must be duplicated at intervals of W = 2t:R in 
the y direction. We choose, without loss of generality, 
^0 ^ Vq ~ 0. Therefore the stresses it generates at {x, y) 
are given by: 



cylinder \ ^ 

XX / J 

k— — oc 



^Abh[y ~ kW][3x'^ + {y- kWf 
[x2 + {y- kWYY 



cr. 



cylinder 
VV 



E 



Abhylx"^ -{y- kWf 
[x2 + {y- A;W)2]2 



cyhnder _ Abhxjx^ ^ {v - kWf 



'xy 



k— — oo 



[x^ + {y- kwy 



(4) 



(5) 



(6) 



/(^) 



-Abh 



{y/W + z)[^x/WY + {y/W + zf 
W[{x/Wf + {y/W + zff 



(10) 



Since f[z) ^ 1/z at large distances from the origin, we 
can deform the contour so that it captures only the poles 
of f{z) (the integral on the circle at infinity vanishes even 
though the decay is only ^ l/z, due to the cot(7rz) term). 
The function f[z) can be rewritten as: 



f{z) - ~{Abh/W) 



{y/W + z)[?,{x/Wf + {y/W ^ zf 
(z + (y + ix)/WY-{z + {y- ix)/Wf 



(11) 

Clearly, it has two poles, each of order 2. Upon summing 
the residues we find: 



(Jxx = ^ g{z)f{z)dz = 

iAbhu'^x 2, , ■ \ i^^T\ 
-csc^(7r(y - ix)/W) 



2W^ 
-KAbh 

' 2W 



cot(7r(y + ix)/W) + C.C. 



(12) 



In a similar fashion one obtains: 



iAbhn^x 



G. 



IX Abh 



csc^(7r(y - ix)/W) 



2W 



cot(7r(y4^ix)/W) + aC., 



(13) 
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and 



^-v = -^Pf7^^^c'('^(y - + C.C. (14) 

Upon integration, one obtains Eq. (5) in the main text, 
describing the interaction energy of two dislocations on 
a cyhnder as a function of their separation. The energy 
landscape, with and without an additional constant force 
along the circumference, is shown in Fig. 1. 




-1 -0.5 0.5 1 

(a) sffi 




-1 -0.5 0.5 1 



X 

2nR 

(b) 

FIG. 1: (a) Contours of E(x,y) for the interaction energy of 
two dislocations on an infinite cylinder, as given by Eq. (5) in 
the main text. The y axis is along the circumference, and the 
X axis is along the axis of symmetry, (b) Energy contours of 
E{x,y) — Gy with ^^^^ = 1. As explained in the main text, 
a combination of chemical forces and turgor pressure would 
create such a force G on an edge dislocation with a Burgers 
vector along x. This force creates two saddle points denoted 
by S in the potential energy landscape. If thermal activation 
were the dominant nucleation process, the properties of these 
saddle points would control creation of dislocation pairs 0]. 



Numerical simulation of a system of in- 
teracting dislocations 

To study the dynamics of many interacting disloca- 
tions, we numerically simulate a system of 30 active, 
overdamped edge dislocations moving on the surface of 
a cylinder under the influence of a disordered potential 
created by 10,000 inactive dislocations. The parameters 
are chosen according to the estimates of the biological 
parameters, described in the main text. Their Burgers 
vector ±6 is along the x axis, which coincides with the 
cylinder's long axis, and we insert equal numbers of dis- 
locations with positive and negative b. The initial posi- 
tions are random along the circumference, and uniformly 
distributed along the cylinder on the interval [0,d], with 
d « 15006; for both active and inactive dislocations. The 
force G on a dislocation consists of two terms: one arises 
from the Gpressure = bcTxx term due to the turgor pres- 
sure inside the bacteria, the second due to the energy 
change u associated with the biochemical reactions in- 
volved in the elongation process, as discussed in the main 
text. Another contribution to the force arises from the 
elastic interactions between the dislocations, which (for 
an infinitely long cylinder) are given by Eqs. ([7)1 . ([S]) . ((T^ 
and (fT4|) . The number of inactive dislocations, 7Vi„, is 
10,000 for the most biologically relevant case, shown in 
Fig. 3 of the main text. A cutoff Fc on the maximal 
elastic forces exerted by a dislocation is imposed, associ- 
ated with the short distance cutoff, which we take to be 
of order b, the lattice spacing. All simulations are run at 
zero temperature, and solved by direct integration of the 
equations of motion. The total simulation running time 
is denoted by T,. An important parameter choice in the 
simulations involves the mobility tensor fj.ij defined by 
Eq. (4) of the main text. We take this tensor to be diag- 
onal in the {x, y) coordinate system defined in Fig. 1 of 
the main text, with, ci glide mobility jig — f^xx 

and climb 

mobility /ic — fiyy. We approximate the predominantly 
climb motion set up by the elongation enzymes via the 
choice = 10/^g. The small non-zero glide mobility al- 
lows us to model the fluctuations of the motion along the 
cylinder's long axis, which could arise from disorder in 
the peptidoglycan mesh. 

To illustrate further the effects of elastic interactions 
between dislocations, we eliminated all inactive disloca- 
tions and simulated the dynamics of 30 active disloca- 
tions strongly interacting on a cylinder, again with a 
Burgers vector b ~ ±x . The results are shown in Figs. 
2 and 3, for two values of the driving force G, differ- 
ing by a factor of 10. For the stronger driving field the 
elastic interactions are less important, and one obtains 
ordered circumferential motion. For the weaker field, the 
attractions between dislocation of opposite sign of the 
Burgers vector are strong enough such that they occa- 
sionally overcome the driving force and meet, leading to 
their mutual annihilation. A detailed study of the an- 
nihilation process of two interacting dislocations will be 
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given elsewhere 

The dimensionless parameters for Fig 3 in the 
main text and Figs. 2 and 3 are described in Table 1. 




0.5 3 



FIG. 2: Numerical simulation of the dynamics of 30 active 
dislocations, with strong elastic interactions. The disloca- 
tions are driven by a combination of chemical forces and tur- 
gor pressure in the circumferential direction, proportional to 
the sign of the Burgers vector. The temperature is zero and 
^^^^ = 3. The X and y coordinates of the dislocations were 
chosen randomly and uniformly along the axes, and W/L — 1. 
The climb to glide mobility ratio is /J-c/jJ^g ~ 10. The red and 
blue points mark the starting positions of the ±& dislocations, 
and the red and blue lines correspond to their trajectories. 
The black circles mark the end of the cylinder. A snapshot 
of the simulation is taken after t — 0.07^^^. A video showing 
the dynamics is attached. 
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FIG. 3: Numerical simulation of the dynamics of 30 active 
dislocations, for the same conditions as for Fig. [2] but with 
a larger field ^^^^ = 30. A video showing the dynamics is 
attached. 



Param. Fig 3 (main) Fig [2] (supp.) Fig [3] (supp.) 



N,„ 10,000 

fic/t^g 10 10 10 

FJG 1/3 100 10 

d/W 1 1 1 

Tsfic§ 0.7 0.07 0.7 



Table 1: Simulation parameters for Fig 3 in the 
main text and Figs. [2] and [3] of the Supplementary 
Material. 

All relevant parameters for the simulations are given 
in the table using the relevant dimensionless ratios. In all 
simulations the number of active dislocations is 30, and 
the number of inactive ones, Nin, is varied. W = 27ri? 
is the circumference and d is length of the interval along 
the X axis where the dislocations (active or inactive) were 
randomly chosen. All simulations use a ratio d/W = 1, 
which is biologically relevant. indicates the relative 

strength of the driving field G versus the elastic interac- 
tions, whose strength is given by the core energy Ab^h. 
The fraction Hc/ fJ-g gives the mobility ratio between climb 
and glide, as explained above. Fc determines the maxi- 
mal elastic force created by the interactions of two dis- 
locations, when their separation is of the order of the 
lattice spacing b. For Fig 3 in the main text, all parame- 
ters arc chosen to be the biologically relevant ones. Figs 
2 and 3 of the Supplementary Material have unrealisti- 
cally large elastic interactions, and are used to illustrate 
the physics of the elastic interactions and their effects on 
the dynamics. 



Characterizing the quenched disorder 

As discussed in the manuscript, inactive dislocations 
interact with the active ones via long ranged elastic 
forces, thus creating a quenched disordered background. 
To estimate the typical magnitude of the disorder, we 
consider a typical E. coli cell with a length of 3 and 
a circumference of 3 fim, and assume an average glycan 
strand length of 30 units, with a unit cell size of 2 nm. 
This implies approximately 10,000 inactive dislocations. 
To find the typical force associated with them, we con- 
sider 10,000 edge dislocations with a Burgers vector along 
the bacteria's long axis (with an equal number pointing 
in the positive and negative direction), and with a ran- 
domly chosen position on the cylinder. We then calculate 
the axx and (Jxy components of the stress tensor at ran- 
domly chosen points. Fig. 4 shows the distribution of 
the resulting forces F^. and Fy exerted by this stress ten- 
sor on an edge dislocation. The standard deviation of 
these forces are found to be of the order of 10^^ dyne for 
both axx and a^y ■ Interestingly, this is the same order of 
magnitude as the force due to the turgor pressure and the 
biochemical forces, estimated in the main text. Hence, 
we expect that the interactions with the background of 
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F (dyne) xio ° 

FIG. 4: Numerical determination of the force distribution 
due to the stress created by = 10, 000 inactive dislocations 
randomly positioned on a cylinder of length L = 3fim and 
with a circumference of Vl^ = 3^m. All inactive dislocations 
were taken as edge dislocation with a Burgers vector along the 
bacterium's long axis, with an equal number of dislocations 
for each sign of the Burgers vector. 100 disorder realizations 
were generated, and for each the total stress was evaluated 
for 10,000 randomly chosen points, using Eqs. p2|) and (|14[) . 
with Ab^h = 2 ■ 10~13 erg. The figure shows the distribution 
of the X and y components of the forces exerted by this stress 
tensor on an edge dislocation with b — x, using Eqs. (O and 



inactive dislocations can be significant, and could pin the 
motion of the active dislocations, as illustrated in Fig. 3 
of the main text. 



Estimation of the biological parameters 

We estimate the relevant biological parameters for the 
case of E. coli, to isolate the important physical pro- 
cesses. 

Dislocation core energy: The core energy associated 
with a single dislocation, Ec, arises from broken bonds 
and localized elastic disturbances, effects which we now 
estimate. According to Rcf. the Young modulus of 
E. coli is anisotropic, with the elastic constant differ- 
ing by roughly a factor of 2 between the two symmetry 
axes (of the NAG-NAM chains and of the peptides 0)- 
The 3d Young modulus is of the order y = 2 ■ l{f N/m?, 
with similar values from Ref. 0], which also discusses 
non-linear effects Q . Upon assuming the sacculus is ap- 
proximately incompressible, we have a 3d Poisson ratio 
V = 1/2, and therefore Y = 2/i(l + v) = 3/i. In this limit 
we find A ^ in Eq. (3) in the main text. The Burg- 
ers vector h and the cell wall thickness h should both 
be of order several nm for gram negative bacteria like 
E. coli. With b h 2nm, we find that the contri- 
bution of the elastic forces to the core energy is of the 



create a dislocation, we must also take into account the 
energy of breaking the peptide bonds. Upon estimating 
this energy as 3 kcal/mole, we find an additional con- 
tribution of ^jpept^de _ 2 • 10-"erg. Thus, we expect 
Ec ^ 4:- 10~^^erg. This value is significantly higher than 
ksT, which is approximately 4 • 10~^'^erg at 37 degrees 
Celsius. 



Interactions between dislocations: The force due 
to a dislocation scales as hAb^ ■ 1/ distance = 2 ■ 
erg / distance. Let us find at which distance the 
interaction force between two dislocations would be com- 
parable to that arising from the osmotic pressure. The 
latter seems to be of the order of an atmosphere in 
gram-negative bacteria Q , which implies a stress tensor 
(T ^ pR ^ lO^Pa, and a Peach-Kohler force h-b-cr ^ 10~® 
dyne. The resulting distance is of the order of a nanon- 
meter, suggesting that G > G* as discussed above. If 
the main force on a dislocation arises from the turgor 
pressure, the Peach-Kohler force may help overcome in- 
teractions, unless cells are starved, enzymatic machinery 
is disabled, etc. In healthy cells, with active glycan exten- 
sion machinery, turgor pressure alone might be sufficient 
to drive the unbinding and subsequent motion of dislo- 
cations. There is, however, an additional complication: 
experiments suggest that the density of inactive disloca- 
tions far exceeds that of active dislocations, nac <S 
since measured strand lengths are typically tens of repeat 
units [13, Hi- This allows us to estimate the number 
of inactive dislocations as NinacUve ~ 10, 000, where we 
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took an average strand length of 30, and a typical separa- 
tion between adjacent strands of 2 nm. As shown in the 
section "Characterizing the disorder", for this scenario 
the typical net force due to this quenched potential is of 
the order of 10~^ dyne, comparable to the "driving force" 
G. Thus, we expect elastic interactions with the inactive 
dislocations can be important. This effect is illustrated 
in Fig. 2 of the main text, showing that in this regime 
the interaction with the inactive dislocations can pin the 
motion of the active ones. For weaker driving forces, we 
expect that this effect will be even more pronounced, and 
could significantly reduce the growth rate. The relatively 
short glycan strand lengths measured in experiments can 
be reconciled with relatively long measured proccssivitics 
in two way: Long glycan strands can be cut by enzymes 
after they have been inserted by the elongation machin- 
ery, as described by the term 74 in the rate equations. 
This could be potentially useful to the cell wall growth 
process, as a source of inactive dislocations which can 
later be rendered active. Alternatively, it could be that 
the elongation machinery occasionally 'misses' an inser- 
tion of a NAM or NAG unit, thus creating a break in the 
glycan strand, without having the elongation machinery 
detach from the cell wall. 
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